Instability of two interacting, quasi-monochromatic waves in shallow water 
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We study the nonlinear interactions of waves with a doubled-peaked power spectrum in shallow 
water. The starting point is the prototypical equation for nonlinear uni-directional waves in shallow 
water, i.e. the Korteweg de Vries equation. Using a multiple-scale technique two defocusing coupled 
Nonlinear Schrodinger equations are derived. We show analytically that plane wave solutions of 
such a system can be unstable to small perturbations. This surprising result suggests the existence 
of a new energy exchange mechanism which could influence the behaviour of ocean waves in shallow 
water. 
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The propagation of multiple wave-train systems in 
shallow water has historically received less attention than 
the propagation of one single wave-train. Nevertheless, 
experimental studies carried out by Thompson pj] in rep- 
resentative sites near the coasts of the United States re- 
veal that in 65% of the analyzed data, ocean wave spectra 
show two or more separated peaks in the frequency do- 
main [0, Qj. In this framework, experimental work in 
the laboratory has been performed by Smith and Vin- 
cent [51 . They propagated irregular wave trains with two 
distinct spectral peaks in a wave flume with 1:30 slope 
for different values of the peak frequency and significant 
wave height. From a physical point of view, this condition 
mimics the interaction of two wave regimes, a "swell" and 
a "sea" , propagating in the same direction toward shore 
in shallow water. Their major observation was a decay of 
the higher frequency peak along the flume. They hypoth- 
esized three possible explanations for such experimental 
results: (i) resonant interactions among waves, (ii) bot- 
tom friction that acts differently for the two dominant 
wave peaks, (iii) breaking of the shorter waves enhanced 
by the presence of the longer waves. More recently, using 
a higher order Boussinesq model, Chen et al. y] have 
shown that inclusion of nonlinear interactions, without 
invoking bottom friction or wave breaking, is sufficient 
to account for the decay of the high frequency peak . 
Even though these numerical simulations of the Boussi- 
nesq equation qualitatively reproduce the experimental 
results, the basic physical mechanisms of interaction of 
wave trains with double peaked spectra in shallow water 
are far from being completely understood. 

In this paper we discuss a fundamental instability 
that occurs between two quasi-monochromatic interact- 
ing wave-trains. Furthermore, the focus herein is not to 
attempt to model ocean waves but instead to study lead- 
ing order effects using the simplest weakly nonlinear and 
dispersive model in shallow water, i.e. the Korteveg-de 
Vries (KdV) equation. A great deal of progress in un- 
derstanding wave propagation in shallow water has been 
achieved by investigating the KdV equation, which can 
be considered as the basic weakly nonlinear model for 
unidirectional shallow water waves. The analytical prop- 



erties of KdV (it is integrable by the Inverse Scatter- 
ing Transform B) have improved basic knowledge of the 
nonlinear dynamics of water waves || , |(| , , § . In par- 
ticular after the seminal work by Zabusky and Kruskal 
JlfH , extensive studies have been carried out on the evolu- 
tion of a sine wave in shallow water (see for example [0 ) . 
However, less attention has been paid to the evolution of 
two monochromatic waves with different wave-numbers. 
The subject of this Communication is therefore an in- 
vestigation of the basic nonlinear interaction that may 
take place between two separated narrow-banded wave 
spectra in shallow water. To this aim, under suitable 
assumptions, using a multiple-scale technique, we derive 
from the KdV equation a system of two Coupled Nonlin- 
ear Schrodinger (CNLS) equations. Plane wave solutions 
of such a system are then studied analytically by a stan- 
dard linear stability analysis, resulting in the presence of 
an instability region. 

The KdV equation can be formally derived from the 
Euler equations for water waves @,§ under the assump- 
tion that waves are small (but finite amplitude) and long 
when compared with the water depth at rest. In a frame 
of reference moving with the velocity cq — \fgh, where h 
is the water depth and g is gravity acceleration, the KdV 
equation in nondimcnsional form reads: 



rj t + [ir\r\ x + \rj xxx = 0. 



(1) 



Here r/ = rj(x, t) is the free surface elevation, x and t 
are space and time variables; fj, and A are the nonlin- 
ear and dispersive small parameters: fi — 3a/2h and 
A = (h/l) 2 /6, with a a characteristic wave amplitude 
and I a characteristic wavelength. We are interested in 
investigating the interaction of two waves, centered at 
nondimensional wave- numbers k\ and k-i , propagating in 
the positive x direction. We consider the case of narrow- 
banded spectra, i.e. Aki/ki <C 1, with i = 1,2, Afc^ 
being the characteristic width of spectra around each 
peak. The approach used in our derivation of the CNLS 
equations resembles the one used by R. Grimshaw et al. 
flllf to derive a single higher order NLS equation starting 
from an extended KdV (cubic nonlinearity is included). 
We introduce the following slow space and time variables 
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X = ex and T = et, with e«l, and perform a formal 
expansion of r\ — r](x, t, X, T): 

r)(x,t,X,T) = ^[A(X,T)e mi + B(X,T)e l&2 + c.c] 

+ ^-[A 2 (X,T)e 2iei + B 2 (X,T)e 2l&2 

+ C(X, T)e j(01+e2) + D(X, T)e^ @1 - @ ^ + c.c] 
+ e 2 M(X,T) + (2) 

where 0j = kix — w{t, M(X,T) is a real quantity to 
be considered as a mean flow and c.c. indicates complex 
conjugate. From a physical point of view this represen- 
tation corresponds to a double expansion around wave 
numbers k\ and k 2 . Note that terms at order e 2 include 
the complex envelopes A 2 and B 2 for the higher harmon- 
ics and the complex envelopes C and D for the harmon- 
ics given by the sum and difference of the fundamental 
wave-numbers . 

After substituting the expansion (Q) into Equation (Q) 
and collecting terms for different harmonics, after some 
lengthy but straightforward algebra, we obtain a set of 
equations for the complex envelopes and the mean flow 
M . At order e, the equations for A or B provide the 
linear dispersion relation: = — Xkf (i = 1,2). For 
the second harmonics, at order e 2 , the following relations 
hold between complex envelopes: 



A 2 
C 
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At order e 3 , the mean flow M is related to the envelopes 
A and B as follows: 



Mt + t( \A\ 



\B\ X 



0. 



(4) 



Using Equation ([|) and neglecting terms of order higher 
than e 3 , the equations for A and B then read: 



e 2 (A T -3\klA x )+ie 3 



(5) 



i 2 1 



= 0, 



e 2 (B T -3Xk 2 B x )+ie 3 

fik 2 M 



3Xk 2 B X x + 
2 1 



(6) 



24A k 2 1 



\B\ 2 IB 



0. 



Note that the arguments used to derive Equations (^Hg) 
are not valid when the two carrier waves k\ and k 2 are 
equal (additional care in the expansion should be taken 
to study the case of ki —k 2 ). Considering that at e 2 the 
following relations hold, 
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tBt, 



(7) 



the mean flow can be directly related to the complex 
envelopes A and B: integrating Equation (^) in time we 
get 



M = - 



A 4 (\A\ 2 , \B\ 



12AV k\ k 2 



(8) 



After substituting Equation (0) in (0), (^|), and rescal- 
ing the variable as follows: T — et/3X, X = ex, 
A' = Aefci/i/(A6v^), B' = Bek 2 n/{X&^/2), two defo- 
cusing Coupled Nonlinear Schrodinger equations are ob- 
tained: 



• 1/ -/.'f-'l, +■ '/m. 1,., -i T ^\A\ 2 A-2i'^\B\ 2 A = 0, (9) 



k 2 



B t - k 2 B x + ik 2 B xx ~ iT^\B\ 2 B - 2i^\A\ 2 B = 0, 

K 2 K 1 

(10) 

where we have dropped the primes for simplicity. In 
these new variables the complex envelopes appearing 
in (PHlOj) are proportional to the Ursell number J5|,||(J 
(U r = kaj (kh) 3 ), therefore to the nonlinearity of the sys- 
tem. If we set B = 0, the system reduces identically to 
the well known defocusing Nonlinear Schroedinger equa- 
tion for the complex envelope A, originally derived from 
the Euler equations in the limit of shallow water, for 
narrow-banded spectra and small steepness The 
single defocusing NLS equation has also been obtained 
from the KdV equation in lll| , see also jlj . For a single 
wave train the expansion (0) is nothing but the second- 
order modulated Stokes wave, i.e. the Stokes series in 
shallow water. It is well known that the plane wave so- 
lution of the NLS equation in shallow water (defocusing 
NLS) is stable to side-band perturbations. This contrasts 
to the Benjamin-Feir instability of the focusing NLS that 
describes the propagation of waves in deep water. 

Systems of equations similar to (^-10) are not new in 
various fields of physics. For example, focusing CNLS 
equations have been derived in |l2] ], |l| for optical fibers. 
For particular values of the coefficients of the nonlinear 
terms in (p[]l0|) , the integrable Manakov system is recov- 
ered [[li) (see also jjl§). Concerning water waves, focus- 
ing CNLS equations have been derived for the first time 
by Roskes Q in infinite water depth. CNLS systems de- 
scribing the interaction of two counterpropagating waves 
are discussed in p0[ , Q. Unlike the work by Roskes, the 
attention here is focused on the shallow water limit, i.e. 
in the case where the system reduces to two defocusing 
NLS equations. 

We now consider the linear stability analysis of the 
system (|9|-^0|). Consider a plane wave solution of CNLS 
of the form: 



A = Ae 



B = Be~ 



with 



u JA ^-{\A\ 2 /kl + 2\B\ 2 k 1 /kt), 
u B =-(\B\ 2 /kl + 2\A\ 2 k 2 /ki), 



(11) 

(12) 
(13) 
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i.e. the nonlinear frequency correction in the Stokes ex- 
pansion. The solutions (|ll| ) are then perturbed as fol- 
lows: 

A = A(l + a)e- t{uAt+a \ 
B = A{l + b)er l{uBt+f) \ (14) 

where a, b, a and (3 are small perturbations in amplitude 
and phase. We substitute relations ( [l4| ) into the CNLS 
equations and obtain a system for a, b, a and j3\ after 
linearization, the standard Fourier technique is used for 
solving it: 

a = he l ^ Kx - nt \ b = be li - Kx ' nt \ 
a = ae l{Kx - nt) , (3 = p e l{Kx - nt) , (15) 

where a, 6, a and (3 are constant, 17 = Q(K) is the disper- 
sion relation for perturbed wave- number K . After some 
algebra, the dispersion relation results in a fourth-order 
polynomial function in the variable f2: 

k 2 k\ {ik\K + n) 2 - k 2 K A - 2\A\ 2 K 2 /k 2 ^j x (16) 

({k 2 2 K + n) 2 - k 2 K* - 2\B\ 2 K 2 /k^j = 

16\A\ 2 \B\ 2 K 4 . 

Equation (^|) provides the dispersion relation for the 
perturbation: complex roots, which are functions of the 
parameters ki, ki, A, B and K, originate instability. 
Note that if we set one of the two amplitudes equal 
to zero (for example B = 0), the right hand-side term 
in ([if]) vanishes and roots can be found easily. In this 
case the roots are always real: a single monochromatic 
wave in shallow water is stable to side-band perturba- 
tions. Equation ([!(]) has been investigated by applying 
the exact formula for the solution of a fourth-order poly- 
nomial equation. Complex roots are then found. Results 
are presented in the following way: we fix the nondimen- 
sional wave-number k\ — 1 (this choice is natural if we 
take the characteristic length I to be l/k, with k the 
dimensional wave-number of the first peak in the spec- 
trum) ; amplitudes A and B are then selected in order to 
fix the nonlinearity of the initial conditions. We then 
compute the largest imaginary solution of the polynomial 
function Q and plot it as a surface for different values of 
ki and perturbation K; as ki is increased, the amplitude 
B is kept constant. In dimensional variables this would 
correspond to increasing the wave-number ki and at the 
same time maintaining constant its Ursell number. 

In Fig. [l] we show the instability region for A = 0.2, 
B = 0.2. The plot clearly exhibits an unstable region 
with a growth-rate different from zero. Note that on the 
horizontal axes ki starts from values slightly larger than 
1, because, as has been previously stated, the CNLS sys- 
tem that we have derived is not valid when ki ~ k\. In 
Fig. H we show the same diagram as in Fig. |l|, with 
larger values of the amplitudes A = 0.4 and B = 0.4. 




FIG. 1: Growth-rate as a function of nondimensional wave- 
number &2 and perturbation K for A = B — 0.2 and fci = 1. 
Axes are in nondimensional units. 




FIG. 2: Growth-rate as a function of nondimensional wave- 
number &2 and perturbation K for A = B — 0.4 and ki = 1. 
Axes are in nondimensional units. 



The instability and the growth-rate have now both in- 
creased. We recall that the amplitudes A and B have 
been scaled by a combination of the steepness, of the 
nonlinear and dispersive parameters, respectivelly \x and 
A, of the KdV equation. From a physical point of view, 
increasing the amplitude A or B implies an increase of 
the nonlinearity of the wave system that can be achieved 
for example also by a decrease in the water depth. As the 
the water depth decreases the instability region and the 
growth-rate increase and therefore, double-peak spectra 
are more likely to evolve naturally into single-peak spec- 
tra which are stable in shallow water. Results can be 
summarized as follows: in shallow water the dynamics of 
two wave trains can be unstable; as expected the growth- 
rate and the size of the instability region depend on the 
nonlinearity of the system. 

The derived CNLS system represents a very crude sim- 
plification of the real problem: effects related to a non 
constant water depth could be considered; directionality 
could also be simply included by applying the multiple- 



4 



scale method to the Kadomtsev-Petviashvili equation |4[] 
(an extension of the KdV equation that includes the dy- 
namics of transverse perturbations); higher order effects 
could also be investigated. 

We believe that the present work offers new perspec- 
tives for understanding the dynamics of double-peaked 
spectra in shallow water. Accurate comparison with 
experiments and numerical simulations are definitely 
needed in order to validate the obtained theoretical re- 
sults. As a final remark we would like to stress that 
we have started the derivation of the CNLS equations 



from the KdV; the Inverse Scattering theory furnishes 
a unique method for investigating all its solutions. It 
would be therefore interesting to interpret the unstable 
solutions of the CNLS equations in terms of the Inverse 
Scattering modes for the KdV equation. 
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